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Abstract 

In this paper we present numerical and analytical studies on the complete synchronization 
phenomena exhibited by unidirectionally coupled two variant of Murali-Lakshmanan-Chua cir¬ 
cuits. The transition of the coupled system from an unsynchronized state to a state of complete 
synchronization under the influence of the coupling parameter is observed through phase por¬ 
traits obtained numerically and analytically. 
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1 Introduction 

The phenomenon of chaos synchronization in coupled chaotic systems has been studied extensively 
for the past two decades with a motivation to understand the coherent dynamical behaviour of 
coupled systems. Chaos synchronization in nonlinear electronic circuits has been observed in a va¬ 
riety of circuits [?,?,?,?,?,?,?] and finds potential applications in secure communication. Several 
synchronization phenomena such as complete, phase and lag synchronization have been identified 
in identical and non-identical chaotic systems and studied both experimentally and numerically. 
A detailed study on characterizing the above said synchronization phenomena was done numeri¬ 
cally [?]. Complete synchronization being the strongest among the other types of synchronization 
occurs when the states of the coupled systems converge, irrespective of the mismatch in initial 
conditions. Complete synchronization phenomena exhibited by coupled second order dissipative 
electronic circuits was studied experimentally, numerically [?,?] and analytically [?]. The synchro¬ 
nization of coupled Chua circuits and coupled Murali — Lakshmanan — Chua{MLC) circuits was 
studied numerically [?]. The MLC circuit being the simplest diisipative non-autonomous chaotic 
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circuit was introduced by Murali [?]. The period doubling dynamics of the circuit leading to chaotic 
motion was studied both experimentally and numerically [?,?]• An explicit analytical solution to 
the normalized circuit equations of the MLC circuit were presented [?,?,?]. The MLC circuit 
has the picewise linear nonlinear element, the Chua’s diode, which is linear within the three re¬ 
gions. Since the circuit equations describing the MLC circuit in each of the piecewise linear region 
is a second order differential equation, the equation is solved for each of the linear regions and 
matched across the boundaries. Similar solutions were given to some simple second order nonlinear 
chaotic circuits with piecewise nonlinear element as the active circuit component [?,?,?,?] and for 
some circuits exhibiting Strange Nonchaotic Attractors (SNA) in their dynamics [?]. The complete 
synchronization phenomenon exhibited by undirectionally coupled two MLC circuits was studied 
analytically through phase portraits generated from the analytical solution [?]. 

The variant of the MLC circuit {MLCV) is a simple forced parallel LCR circuit with the Chua’s 
diode as the nonlinear element, connected parallel to the capacitor, was suggested by Thamil- 
maran [?]. The quasiperiodic route and the reverse period doubling route to chaos exhibited by 
MLCV circuit has been extensively studied [?,?]. Further an explicit analytical solution to the 
state variables of the normalized circuit equations are obtained. The solutions thus obtained were 
used to explain the dynamics of the circuit through phase portraits of the state variables. The 
complete synchronization phenomena exhibited by MLCV circuit was studied experimentally [?]. 
However the dynamics of the coupled MLCV circuit has not been yet studied numerically and 
analytically. In this paper, we study the complete sychronization phenomenon exhibited by uni- 
directionally coupled two variant of Murali-Lakshmanan-Chua circuits numerically and provide an 
explicit analytical solution for the state variables of the normalized circuit equations characterizing 
the coupled system. 


Circuit Equations 

The variant of the Murali-Lakshmanan-Chua circuit is a simple forced parallel LCR circuit with the 
Chua’s diode as the only nonlinear element, connected parallel to the capacitor C. The normalized 
state equations of the coupled circuit is given as. 


X = fisin{ujit) - X - y - g{x), 
ij = (3x 


The piecewise linear function g{x) representing the Chua’s diode is given by, 


9{x) 


hx -\- {a — b) if X > 1 
< ax if |x| < 1 

bx — {a — b) if X < — 1 


(la) 

(lb) 


( 2 ) 


where, /3 = {C(LC^), a = CJC, b = Cb/C, h = (Fi/3/Hp), = {VLiC/G) and G = l/R. The 

circuit parameters takes the values C = lO.lSnT, L = 445mFf, R = 1475H and the parameters of 
the Chua’s diode are chosen as Ga = —0.76m5, Gb = —0.41mS' and Bp = IV. The frequency of 
the external periodic force ni is fixed as ui = kl(2iT = 1.116kHz. The MLCV circuit explained 
above acting as the drive system is unidirectionally coupled to another M LCV circuit which acts as 
the response system. The drive and the response systems operate with different initial conditions, 
different values of amplitude of external periodic forcing, are coupled by a linear resistor and a 
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Buffer 



Figure 1: Schematic circuit realization of two identical MLCV circuits with unidirectional coupling 


buffer. The buffer acts as a signal driving element which isolates the drive system variables being 
affected by the response system.The schematic diagram of the coupled MLC circuits is as shown 
in Fig.l. The normalized state equations of the response system is given as, 


X = f 2 sin{u} 2 t) — X —y — g{x ) + t{x — x ), 
y = I3x 

and the piecewise linear function g{x') is given by, 


9{x') 


bx + {a — b) 
/ 

< ax 

bx — {a — b) 


if > 1 

if |x^| < 1 
if < — 1 


(3a) 

(3b) 


(4) 


where, /2 = {F2f3/Bp), 1 V 2 = {^ 2 C/G) and e = (R/Rc) is the coupling parameter. The circuit 
parameters of drive and response systems are fixed as C = 10.15nF, L = AAbmH, R = 14750 and 
the parameters of the Chua’s diode are chosen as Ga = —0.76mS, Gf, = —0.41m5 and Bp = IV. 


2 Numerical Analysis 

In this section, we provide the numerical simulation results of the unsynchronized and synchronized 
states of the coupled system. Further, we present the eigenvalues of the difference system, obtained 
from the normalized state equations of the system, as a function of the coupling paramter. 

From the normalized state equations of the response system given by Eqs.(3) we could infer that 
the two circuits are independent of each other when the coupling resistor Rc takes very larger 
values such that the coupling parameter e = R/Rc = 0. Hence, the coupling resistor blocks the 
response being controlled by the drive and makes it independent of the drive system. The frequency 
of the external periodic force are fixed as uji ^2 = 0.105 while their correponding amplitudes are 
fixed to have different values as fi = 0.38 and /2 = 0.381. In addition to that, the drive and 
response systems are operated with different set of initial conditions as xq = — 0.5,yo = 0.1 and 
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Xq = 0.5, i/q = 0.11. It is to be expeceted that the chaotic attractors of the drive and response 
systems must be unsynchronized because of the mismatch in initail conditions and external forcing. 
Fig.2(a), (b) shows the unsynchronized state of the drive and the response systems in the x — x 
phase plane and the trajectory plot m. x — x plane respectively for e = 0. However, for small 
values of Rc, e = R/Rc > 0, the coupling resistor allows signal flow from drive to response and the 
two systems are strongly coupled. Owing to unidirectional coupling, the dynamics of the drive is 
unaltered while that of the response system varies with the coupling parameter e and hence it is 
controlled by the drive. Fig.2(c),(d) shows the synchronized state of the two systems in the x — x 
phase plane and the trajectory plot in x — x plane respectively for e = 1. It is found that the 
drive and the response systems which are unsynchronized for low values of coupling paramter e are 
completely synchronized for higher values. 


1 

0 





Figure 2: Numerically obtained chaotic attractors in (a) x — y (b) x — sin{u}it) phase planes. 


Stability Analysis 

From the numerical results obtained, we observe that the dynamics of the response is influenced by 
the drive. Hence, the eigenvalues of the response system must depend upon the coupling parameter. 
Since the drive and response systems are piecewise linear, each picewise linear regions of the two 
sytems could be coupled together to get a new set of equations. In this section, the dependence of 
the eigenvalues and their stability on the coupling paramete for each of the piecewise linear regions 
is presented. The difference system obtained from Eqs.(I) and (3) are 

X* = fisin{oJit) — f 2 sin{ijJ 2 t) — x* — y* — (y(x) — g{x )) — ex*, (5a) 

y* = I3x* (5b) 

where x*={x — x'), y*={y — y) and g{x) — g{x) = g{x*) takes the values ax* or hx* depending 

upon the corresponding region of operation of the drive and response systems. Since the circuit 

equations are piecewise linear, the equations representing the identical regions of the two circuits 

are coupled to get a new set of equation for the corresponding regions. The state variables of this 

new set of equations, x*{t), y*{t) are given as, 

* ' 

X = X — X 

y* = y- y\ 
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(6a) 

(6b) 






Figure 3: (a) Unsynchronized motion in {x — x) plane for e = 0.0, /i = 0.38 and /2 = 0.381, 
[xo=-0.5, yo=0.1] and [xq = 0.5, yQ=0.11]. (b) Trajectory plot in (x — x ) plane, (c) Synchronized 
motion in (x — x ) plane for e = 1, /i = 0.38 and /2 = 0.381, [xo=-0.5, yo=0.1] and [xg = 0.5, 
2/0=0.11]. (d) Trajectory plot in (x — x ) plane. 


provided the drive and the response systems exists in identical piecewise linear regions. One can 
easily establish that a unique equilibrium point (xq, 2 /q) exists for equation (6) in each of the fol¬ 
lowing three subsets. 


ii;i = {(x*,2/*)|x* >i}p; = (o,o), \ 

D* = {{x*,y*)\\x*\<l}\O* = {0,0), i (7) 

71*1 = {(x*, 2 /*)|x* <-1}|P: = (0,0),J 

Naturally these fixed points can be observed depending upon the initial conditions Xg and Uq of 
Eqs.(5) which are inturn obtained from the initial conditions of the drive and response systems as, 
Xg = xo — Xg and Uq = yo — ?/g. In fact, Eqs.(5) can be integrated explicitly in terms of elementary 
functions in each of the three regions Dq, and the resulting solutions can be matched across 
the boundaries to obtain the full solution. 

The stability of the fixed points given in Eq.(7) can be calculated from the stability matrices. 
Hence, in the first case, g{x) and g{x ) takes the values ax and ax respectively, corresponding 
to the central region in the v — i characteristics of the nonlinear element, which is taken as the 
Dq region of the difference system. In the Dq region the stability determining eigen values are 
calculated from the stability matrix, 

— (a -|- e -|- 1) 
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The eigen values mi and m 2 are found to be a pair of complex conjugates for e < 0.5682135 while 
they are real and distinct for e > 0.5682135. 

In the second case, g{x) and g{x ) takes the values bx ± {a — b) and bx ± (a — b) respectively, 
corresponding to the outer regions in the v — i characteristics of the nonlinear element, which is 
taken as the regions of the difference system. In the regions, the stability determining 
eigen values are calculated from the stability matrix. 


({b + e + 1) -A 

V /? 0; 


(9) 


The eigen values mi and m 2 are found to be a pair of complex conjugates for e < 0.0519135 while 
they are real and distinct for e > 0.0519135. The eigen values of the difference system in both 
the regions are thus determined by the value of the coupling parameter e. The eigenvalues of the 
difference systems in both regiona sre summarized in Table 1. 


Region 

Fixed Point 

e 

Eigenvalues 

Stability 

n* 

(0,0) 

e < 0.5682135 

e > 0.5682135 

Complex conjugates 
Real and distinct 

Stable focus 
Stable node 

D* 

^±1 

(0,0) 

e < 0.0519135 

e > 0.0519135 

Complex conjugates 
Real and distinct 

Stable focus 
Stable node 


Table 1: Stability of the fixed points in Dq and regions 




Figure 4: Eigenvalues as functions of the coupling parameter e in (a) Dq region (b) D'^^ region 


Fig.4(a) & 4(b) shows the eigenvalues of the Jacobian matrix given by Eq.(8) and (9) in the 
Dq and regions respectively as functions of the coupling parameter e. The red and green lines 
show the two real roots while the blue line shows the real part of the complex conjugate roots as 
functions of the coupling parameter. In the next section, we present explicit analytical solutions for 
the dynamics of the response system driven by the drive system for value of the coupling parameter 
(e > 0) leading to complete synchronization of the drive and response. 
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3 Explicit analytical solutions 


In this section, we present explicit analytical solutions for the unsynchronized (e = 0) and synchro¬ 
nized states (e > 0) of the coupled systems. When the coupling paramter e = 0, the two systems 
are independent of each other. Hence the drive and the response systems given by equations (1) 
and (3) have the same solution for their state variables for all the three piecewise linear regions. 

Analytical solution for e = 0 

An explicit analytical solution for the single or uncoupled MLCV circuit is presented and a chaotic 
attractor resembling that of the numerically obtained one as shown in Fig. (2) is observed in the 
X — sin{ujt) phase plane from the analytical solutions thus obtained [?]. The solutions of the 
normalized equations in terms of the state varaibles y{t) and x{t) obtained for each of the three 
piecewise linear regions Dq, D± are summarized as follows. 

Region : Dq 

Since the roots mi ^2 in this region are a pair of complex congugates, the fixed point (0, 0) cor¬ 
responding to the Dq region is an unstable spiral fixed point. The state variables y{t) and x{t) 
are 

y{t) = e^^{Cicosvt + C 2 sinvt) + Eisin{LOit) + E 2 Cos{u!it) (10) 

x{t) = 1(S) (11) 

Region : D±i 

The roots mi ,2 in this region are also a pair of complex conjugates. Hence the fixed points ki = 0, 
/c 2 = ±(a — h) corresponding to the D±i region is a stable spiral fixed point . The state variables 
y{t) and x{t) are 

y{t) = e^^{Czcosvt -\- C^sinut) + E^sin{ujit) + i? 4 Cos(a;it)±A (12) 

= t(!/) (13) 

where A = (b — a) and -|-A, —A corresponds to D+ and D- regions respectively. 

The above solutions of the state variables can be used to obtain chaotic attractors for both the 
drive and response systems. Fig.5 shows the analytically obtained phase portraits for the chaotic 
attractors of the drive and response systems together in the (x — y) and {x — y ) phase planes 
for an amplitude /i = 0.375, /2 = 0.377 and frequency t<;i ,2 = 0.105 of the external periodic force. 
Since the response system is operated with a different set of initial conditions and amplitude of the 
external periodic force, the chaotic attractors of the drive and response systems are non- identical 
and are unsynchronized. 

Analytical solution for e > 0 

When the coupling parameter is increased ie.(e > 0), the response system is controlled by the 
drive and completely synchronizes with the drive for larger values of the coupling parameter. An 
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Figure 5: Analytically obtained chaotic attractors of drive (blue line) and response (red line) 
systems showing the unsynchronized state of the two systems for e = 0. 



explicit analytical solution to the dynamics of the response system could be obtained by finding a 
solution to the normalized state variables of the difference system given by Eq.5. The solution of 
those equations are, [x*{t;tQ,XQ,yQ), y*(t; to) which the initial conditions are written 

as (t, X*,y*) = {to,XQ,yQ). From the solution x*{t) and y*{t) thus obtained the state variables x (t) 
and y (t) can be found using Eqs.(6). Since Eq.(5) is piecewise linear, the solution in each of the 
three regions can be obtained explicitly. 

Region : 

In this region g{x) and g{x') takes the values ax and ax' respectively. Hence the normalized 
equations obtained from Eqs.(5) are, 

X* = —{a + € + l)x*—y* + fisin{uJit) — f2sin(uj2t), (14a) 

y* = fix* (14b) 

Differentiating Eq. (14b) with respect to time and using Eqs. (14a, 14b) in the resultant equation, 
we obtain 

ij* + Ay* + By* = fisin{(jJit) - f 2 sin{uj 2 t) (15) 

where, A = a + e + 1 and B = (3. The roots of the equation (16) is given by 

-{A) ± 

- - 2 -. 


case: a 

Eor e < 0.5682135, (A^ < AB) and hence the roots mi and m 2 are a pair of complex congugates 

given as mi ^2 = u ± iv, with u = and v = xhe general solution to eq.(15) can be 

written as, 

y*{t) = e^^{Cicosvt + C 2 sinvt) + Eisinuit + E 2 Cosu}it + E^sinoj 2 t + EiCosuj 2 t (16) 







where Ci and C 2 are integration constants and 


^ AV + {B-ui^y 

p _ -Aujifif3 

^ AV + iB-coi^y 

^ 712^2^ + {B- u;2^y 

p _ AuJ2f2l3 

* A^U2^ + {B- UJ2^Y 

Differentiating Eq.(16) and using it in Eq.(14b) we get, 

x'(t) = i(r) (17) 

The constants Ci and C 2 in the above equations can be evaluated by solving both Eqs.(16) and 
(17) for Cl and C 2 at a suitable initial instant to, with Xq and Uq as initial conditions at time 
t = to, provided the trajectory of the dynamical system just enters the region Dq at time to- The 
constants Ci and C 2 thus obtained are 


Cl = - {{vyQCOsvto — i/3xQ — ux*Q)sinvto) + {{tdiEi — uE 2 )sinvto — vE 2 Cosvto)cosuJito 

— {{ 1 JJ 1 E 2 + uEi)sinvto + vEicosvto)sinuJito + ((012 Es — uEi)sinvto — vE/icosvto)cosuj 2 to 

— ((a; 2 -E 4 + uEo)sinvto + vE3Cosvto)sinuj2to} 

C 2 = -{((/3xo — uxQ)co-svto + vyQ-sinvto) — {{cdiEi — uE 2 )cosvto + vE 2 sinvto)cosuJito 

+{{ijJiE 2 + uEi)cosvto — vEisinvto)sinuJito — ((012 Es — uE4)cosvto + vE4sinvto)cosu2to 

+(((U 2 -E 4 + uEo)cosvto — vEosinvto)sinu: 2 to} 

From the results of y*{t), x*{t) obtained from eqs.(16), (17) and y{t), x{t) obtained from Eqs.(lO), 
(11) X (t) and y (t) can be obtained from Eqs.( 6 ). 

case: b 

For e > 0.5682135, > AB) and the roots mi and m 2 are real and distinct. The general solution 

to eq.(16) can be written as, 

y*{t) = + C' 2 e™'^* + Eisin{ixit) + E 2 Cos{uJit) + Essin{uj 2 t) + EiCos{iX 2 t) (18) 

where Ci and C 2 are the integration constants and the constants Ei, E 2 , E^, E 4 are the same as in 
case(a). Differentiating Eq.(18) and using it in Eq.(14b) we get, 

x*{t) = ^(y*) (19) 

The constants Ci and C 2 in the above equations can be evaluated by solving both Eqs.(18) and 
(19) for Cl and C 2 at a suitable initial instant to, with Xq and y^ as initial conditions at time 
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t = to, provided the trajectory of the dynamical system just enters the region Dq at time to- The 
constants Ci and C 2 are, 


Cl 


C 2 


g-m2to 

m 2 — mi 


{{/3xq - miyo) + - (jJiEi)cosuJito + (a;i£'2 + miEi)sinuJito 

+{m2Ei — ijj2E‘i)cosuj2to + (W 2 -E 4 + m2E‘i)sinijJ2to} 


g-mito 

mi — m2 


{Oo 


m 2 yl) + (m 2 -E 2 - (jJiEi)cosoJito + (a;i £'2 + m 2 Ei)sinuiito 


+{miEi - oj2E‘i)cosoj2to + (^2-^4 + mi£^3)sina;2to} 

From the results of y*{t), x*{t) obtained from eqs.(18), (19) and y{t), x{t) obtained from Eqs.(lO), 
(11) X (t) and y (t) can be odtained from Eqs.( 6 ). 


Region : 

In this region g{x) and g{x ) takes the values bx±{a — b) and bx ± (a — b) respectively. Hence the 
normalized equations obtained from eqs. (5) are, 

x* = —{b + e + l)x* — y* + fisin{uJit) — f 2 sin{L 02 t), ( 20 a) 

y* = j5x* (20b) 

Differentiating Eq.(20b) with respect to time and using Eqs.(20a,20b) in the resultant equation, we 
obtain 

ij* + Cy* + Dy* = fisin{ujit) - f 2 sin{uj 2 t) (21) 

where, C = 6 + e + 1 and D = j3. The roots of the equation (21) is given by 

_ -(C) ± - ^D) 


case: a 

Eor e < 0.0519135, (C^ < AD) and hence the roots m 3 and m^ are a pair of complex congugates 

given as m 3,4 = uEiv, with u = ^ and v = ^ ^ . The general solution to eq.(16) can be 

written as, 

y*{t) = e^^{Cocosvt + C^sinvt) + E^sinuit + Eocosuit + Ejsinu} 2 t + E^cosijj 2 t (22) 

The constants E^, Eq, Ej, Es are the same as the contants Ei, E 2 , E^, E 4 in case ; a of Dq region 
except that the constants ASzB are replaced with C^D respectively. Differentiating Eq.(22) and 
using it in Eq.(20b) we get, 

X'(t) = i(r) (23) 

The constants C 3 and C 4 are the same as Ci and C 2 in case : a of Dq region except that the 
constants Ei, E 2 , E^, E 4 are replaced with the constants E^, Eq, Ej, Eg respectively. Erom the 
results of y*{t), x*{t) obtained from eqs.(22), (23) and y{t), x{t) obtained from Eqs. (12), (13) 
X (t) and y (t) can be obtained from Eq.( 6 ). 
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case: b 


For e > 0.0519135, (C^ > 411) and the roots m 3 and m 4 are real and distinct. The general solution 
to eq.(24) can be written as, 

y*{t) = C'se’”®* + ( 746 ™"** + E^sin^LOit) + EQCOs{uJit) + Ejsin{uj 2 t) + E 3 Cos{uj 2 t) (24) 

where C3 and C4 are the integration constants and the constants E3, Eq, E^, Eg are the same as in 
case(a). Differentiating Eq.(24) and using it in Eq.(20b) we get, 

x*{t) = ^{y*) (25) 

The constants C3 and C4 are the same as Ci and C2 in case : a of Dq region except that the con¬ 
stants El, E2, E3, E4 are replaced with the constants E5, Eq, Ej, Eg respectively. Erom the results 
of y*{t), x*{t) obtained from eqs.(24), (25) and y{t), x{t) obtained from Eqs. (12), (13) x (t) and 
y (t) can be obtained from Eq.( 6 ). 

Now let us briefly explain how the solution can be generated in the (x* —y*) phase space. Thus 
if we start with the initial conditions x*{t = 0 ) = XQ,y*{t = 0) = yj) in the region Dq region at 
time t = 0, the arbitrary constants Ci and C 2 get fixed. Thus x*{t) evolves as given in Eq.(18) 
or Eq.(21), depending upon the value of e, up to either t = Ti, when x*(Ti) = 1 and x*(Ti) > 0 
or t = T[, when x*{T') = — 1 and x*{T[) < 0. Knowing whether Ti < T[ or Ti > T[ we can 
determine the next region of interest (D^^) and the arbitrary constants of the solutions of that 
region can be fixed by matching the solutions. The procedure can be continued for each successive 
crossing. In this way, the explicit solutions can be obtained in each of the regions Dg, 
However, it is clear that sensitive dependence on initial conditions is introduced in each of these 
crossings at appropriate parameter regimes during the inverse procedure of finding Ti, T^, T 2 , T 2 ,..., 
etc. from the solutions. The analytical solutions obtained above for the response system can 
be used to explain the phenomena of complete synchronization through phase portraits. The 
amplitude and frequency of the drive and the response sytems are fixed as /i = 0.375 ,/2 = 0.377 
and ti;i ^2 = 0.105 respectively. Also the systems are operated with different set of initial conditions 
given by xq = —0.5,yo = 0-1 and Xg = 0.5, yg = 0.11 . Eigure ( 6 ) shows the non-identical chaotic 
attractors of the drive and the response systems with the blue line giving the chaotic attarctor of 
the drive and the red line of the response. Eig .6 shows the Unsynchronized and Synchronized states 
of the coupled systems for two different values of the coupling parameter e. Eigure 6 (a) shows the 
unsynchronized state of the coupled systems for the coupling parameter e = 0 in the (x — x*) phase 
plane and the corresponding trajectory in the x* = x — x plane. As the value of the coupling 
parameter is increased, the response system completely synchronizes with the drive. Eigure 6 (c) 
shows the complete synchronization of the coupled systems in the (x —x*) phase plane for the value 
of the coupling parameter e = 1.0 and the corresponding trajectory in the x* = x — x plane as 
in Eig. 6 (d). Erom the phase portraits and the time series plots obtained it could be inferred that 
for the coupling parameter taking the value e = 1.0, the response system which is operating with 
a different set of initial condition and for a different value of external periodic force, completely 
synchronizes with the drive. 
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Figure 6: (a) Unsynchronized motion in {x — x) plane for e = 0.0, fi = 0.375, /2 = 0.377, [a:o=-0.5, 
yo=0.1] and [xq = 0.5, ?/q=0.11]. (b) Trajectory plot in (x — x ) plane for e = 0. (c) Synchronized 
motion in (x — x ) plane for e = 1.0, /i = 0.375, /2 = 0.377, [xo=-0.5, yo=0.1] and [xq = 0.5, 
2 / 0 = 0 . 11 ]. (d) Trajectory plot in (x — x ) plane for e = 1.0 


4 Conclusion 

In this paper, we presented numerical and explicit analytical studies on unidirectionally coupled 
two variant of MLC circuits exhibiting complete synchronization in their dynamics. The complete 
synchronization of the response system which possesses a non-identical choatic attractor with that 
of the drive, for a particular value of the coupling parameter is studied through phase portraits 
obtained numerically and analytically. 


Acknowledgement 

This work is supported by Alagappa Chettiar College of Engineering & Technology under Technical 
Education Quality Improvement Programme(TEQIP)-II scheme. 


12 










